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Abstract 

A gram-negative aliphatic hydrocarbon-degrading bacterium SJTD-1 isolated from oil-contaminated soil was identified as 
Pseudomonas aeruginosa by comparative analyses of the 16S rRNA sequence, phenotype, and physiological features. SJTD-1 
could efficiently mineralize medium- and long-chain n-alkanes (C 12 -C 30 ) as its sole carbon source within seven days, showing 
the most optimal growth on n-hexadecane, followed by n-octadecane, and n-eicosane. In 36 h, 500 mg/L of tetradecane, 
hexadecane, and octadecane were transformed completely; and 2 g/L n-hexadecane was degraded to undetectable levels 
within 72 h. Two putative alkane-degrading genes (gene 3623 and gene 4712) were characterized and our results indicated 
that their gene products were rate-limiting enzymes involved in the synergetic catabolism of C 12 -C 16 alkanes. On the basis 
of bioinformatics and transcriptional analysis, two P450 monooxygenases, along with a putative AlmA-like oxygenase, were 
examined. Genetically defective mutants lacking the characteristic alkane hydroxylase failed to degrade n-octadecane, 
thereby suggesting a different catalytic mechanism for the microbial transformation of alkanes with chain lengths over C q8 . 
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Introduction 

Oil pollution poses a severe threat to ecosystems, because of its 
refractory resistance to environmental degradation. Many remedia- 
tion technologies have been used for the removal of contaminating oil 
residues [1]. Bioremediation is considered to be a safe, efficient, and 
economically viable alternative to physicochemical methods for the 
elimination of oil contaminants. In the past 30 years, several research 
groups around the world have tried to build robust biocatalysts for the 
efficient removal of contaminating oil residues. Although these 
attempts have not succeeded so far, the constant search for an 
effective biocatalyst has resulted in the discovery of a large number of 
microorganisms that could be used for bioremediation [2,3]. 

Alkanes are the most abundant hydrocarbons (with an estimated 
abundance of 20-50%) in crude oil [4]. As alkanes are non-polar 
molecules with very low chemical activity, their utilization by 
microorganisms faces significant challenges, owing to known 
factors, such as low water solubility, high degree of accumulation 
in cell membranes, and higher activation energies [5,6]. None- 
theless, at least sixty genera of aerobic bacteria, and five genera of 
anaerobic bacteria [7], such as Pseudomonas [8], Acinetobacter 
[9], Rhodococcus [10] and Dietzia [11], have been reported to 
possess the ability to degrade aliphatic hydrocarbons. Among 
them, Pseudomonas was found both in soil as well as in aqueous 
environments [12]. Several Pseudomonas strains are known to 
utilize aliphatic hydrocarbons as its sole carbon sources [13—15]. 
P. aeruginosa RR1 and P. fluorescens CHA0 degrade n-alkanes 
ranging from C 12 to C 34 [16,17]; and P. aeruginosa DQ8 can 



grow in the presence of w-tetrodecane, n-docosane, w-triacontane, 
and w-tetrocontane [8] . 

Although many strains have been reported to utilize hydrocar- 
bons, most of them can only use a narrow range of substrates. For 
example, Bacillus stearothermophilus can only grow in media 
containing hydrocarbons of chain length C 15 to Ci 7 [18], whereas 
A. borkumensis API, SK2, and SK7 can only utilize alkanes 
ranging from C 6 to C 16 [19]. Indeed, very few strains such as 
Acinetobacter baylyi ADP1 and Thermus sp. C2 can degrade a 
wide range of hydrocarbons. However, their genetic characteristics 
remain elusive and not much is known about the mechanisms 
through which these microorganisms break down long chain 
alkanes present in refractory oil residues [20,21]. 

In this study, a P. aeruginosa strain SJTD-1 was isolated from 
oil-contaminated soil; its hydrocarbon utilization capability and n- 
alkane breakdown efficiency were investigated. Based on its whole- 
genome DNA sequence [22], several alkane hydroxylases were 
characterized. These included two AlkB monooxygenases, two 
P450 monooxygenases, and one AlmA-like monooxygenase. 
Varying transcriptional expression levels of these genes induced 
by Ci 2 -C 2 4 alkanes indicated the presence of a complex 
hydrocarbon breakdown mechanism in P. aeruginosa SJTD — 1 . 

Materials and Methods 

Chemicals and media 

n-Decane (>99% pure) was purchased from Alfa Aesar Organic 
Co., Inc (Tianjing); n-dodecane, n-tetradecane, K-hexadecane, 
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and n-octadecane (all >99% pure) from Sangon (Shanghai, 
China); and n-pentadecane, n-eicosane, n-docosane, n-tetraco- 
sane, n-triacontane, and n-hexane (of HPLC gradient grade) from 
Sigma-Aldrich (St. Louis, MO, USA). All other reagents used in 
this study were of analytical reagent grade. 

Luria-Bertani (LB) medium (tryptone 10.0 g, yeast extract 5.0 g, 
NaCl 10 g/L) and basal salt medium (BSM) (K 2 HP0 4 3.815 g, 
KH 2 P0 4 0.5 g, (NH 4 ) 2 HP0 4 0.825 g, KNO s 1.2625 g, Na 2 S0 4 
0.2 g, CaCl 2 0.02 g, FeCl 3 0.002 g, and MgCl 2 0.02 g/L) were 
used in this study. To examine the utilization of alkanes by strain 
SJTD-1, both liquid and solid alkanes, maintained at room 
temperature, were first dissolved in n-hexane to form 500 mg/ mL 
alkane-hexane solutions. These solutions were then added into the 
BSM medium to attain various concentrations. The n-hexane was 
neither toxic to the strain, nor was it utilized by the strain. 

Strain Isolation 

Samples used for bacterial enrichment were collected from the 
oil-contaminated soil in the Daqing Oil Field of China (no specific 
permissions were required and this work did not involve any 
endangered or protected species). Approximately 5 g of the soil 
sample was inoculated into a 250 mL flask with 100 mL BSM 
liquid medium containing 2 g/L hexadecane, and was shaken at 
180 rpm for seven days, at 37°C. A 5 mL culture was then 
inoculated into 100 mL of fresh BSM liquid medium with n- 
hexadecane, and cultured as above. Enrichment cultures after 
rounds of enrichment were diluted and plated onto BSM agar 
plates pre-coated with hexadecane. Bacterial colonies grown on 
plates with different morphology were tested for their n-alkane 
utilizing capabilities. One strain, which exhibited a fast growth 
rate was purified and designated as SJTD-1. 

16S rRNA analysis and phylogenic tree construction 

The newly isolated strain SJTD-1 was identified based on its 
morphological, physiological, and biochemical properties listed in 
Bergey's Manual of Determinative Bacteriology [23]. After its 
genomic DNA was extracted with standard molecular biology 
techniques [24], the 16S rRNA gene of strain SJTD-1 was 
amplified using Bacterial 16S rDNA Kit (TaKaRa Biotechnology 
Co., Ltd. Dalian, China). The PCR program is first denaturized at 
94 °C for 5 min, then denaturized at 94 °C for 1 min, annealed at 
55 °C for 1 min and elongated at 72 °C for 1.5 min; Repeated this 
process for 30 times, and then further elongated at 72°C for 
another 5 min. After the PCR reaction, the fragments were 
sequenced using the primer pairs 1 6S-seq-F/R (Table S 1), and this 
16S rRNA sequence was deposited into GenBank (Accession No. 
JQ95 1926.1). A phylogenic tree based on the 16S rRNA 
sequences of the SJTD-1 strain and some other strains was 
analyzed with MEGA 5.0, using the Neighbor Joining method 
with 1,000 replications. The genetic distances were calculated with 
the Kimura two-parameter distance model. 

Growth rate of SJTD-1 in pure n-alkanes 

The cell growth rate of strain SJTD-1 in n-alkanes was 
determined by the Automatic Growth Curve Analyzer (BioScreen 
Testing Service, Inc., CA, USA). Briefly, a single colony of strain 
SJTD-1 was inoculated into a 100 mL baffled flask containing 
10 mL LB broth, and cultured on a rotary shaker overnight at 
37°C. The cells were then harvested by centrifugation at 
8,000 rpm for 5 min and washed thrice with sterilized water to 
exclude the residual organic compounds. The cells were then re- 
suspended in the BSM medium to OD 6() o~2.0, to form the 
inoculum. Subsequently, 10 |iL of the cell pellets were inoculated 
into the 500 uL well of a 10x10 multi-well plate, which contained 



190 uL of the BSM medium amended with pure n-alkanes (Cio— 
C 24 ) of different concentrations. The initial cell concentration in 
each well at OD 60 o was 0.1. All n-alkanes were prepared as n- 
hexane solutions. Wells containing cells without n-alkanes and 
those containing n-alkanes without inoculum were used as blanks. 
The 10x10 multi-well plates were loaded onto the Automatic 
Growth Curve Analyzer and incubated at 30°C with constant 
shaking (180 rpm) for seven days. The cell densities were 
determined by detecting the OD 600 every hour. All experiments 
were repeated thrice, and the results shown were the average 
values of three replicates, along with the standard errors. 

Degradation analysis of SJTD-1 in the case of pure 
n-alkanes 

The degradation efficiencies of strain SJTD-1 for n-alkanes of 
different carbon lengths (C ]4 -C 16 ) were determined according to 
the loss of substrate from each liquid culture. A SJTD-1 cell 
inoculum was prepared as described above. The cells were 
distributed in flasks containing 100 mL BSM medium, amended 
with 500 mg/L of pure C 10 -C 24 n-alkanes (which was the sole 
carbon source), and cultured for seven days at 30°C in the shaker. 
The initial OD 60 o was 0. 1 . Flasks without cell inoculum were used 
as blanks to assess the abiotic loss. The cultures were taken out at 
different time points, to estimate the cell concentrations and the 
alkane residues. For the residual analysis of alkanes, 3 mL cell 
cultures were collected into 10 mL-glass centrifuge tubes and 
extracted with 1 mL n-hexane for the first time. The mixtures 
were shaken vigorously, centrifuged at 5,000 rpm for 10 min. 
Then, the organic layer was collected and the aqueous layer was 
extracted twice, with 1 mL n-hexane each time. The organic 
extracts were pooled together to a final volume of 3 mL, and dried 
with anhydrous sodium sulfate. For each sample, n-pentadecane 
(100 mg/L) was added before the extraction and used as an 
internal standard. All extractions were performed in triplicate, and 
the results were expressed as average values with standard errors. 

The degradation efficiencies of strain SJTD-1 for n-hexadecane 
of various concentrations (250, 500, 1,000, 2,000 mg/L) were also 
analyzed as above. After strain SJTD-1 was cultured for one to 
seven days, the residual n-hexadecane at different time points were 
extracted and determined with three replicates. 

Analytical methods of n-alkanes concentration 

The concentrations of n-alkanes were determined with the gas 
chromatography-mass spectrometry (GC-MS) technique, using a 
GC/MS system (7890A GC/5975C MS, Agilent, USA) equipped 
with a fused DB-5 MS capillary column (0.25 mmx30 mx0.25 \lm). 
The GC program was listed as below. Helium was used as the carrier 
gas with a constant flow rate of 1 .0 mL/min. The split ratio was 10:1, 
and the temperatures of the injector and the connector were 270°C 
and 280°C, respectively. The temperatures of the ion source and the 
quadrupole were 230°C and 150°C, respectively. The column oven 
temperature was kept at 150°C for 2 min, and then raised to 200°C 
at a rate of 5°C/min, followed by an increase to 290°C, at a rate of 
30°C/min, and at an isotherm of 290°C. The ionization mode was 
set as EI+, 70 ev, and the voltage of the detector was 1,388 V. The 
GC-MS spectra were analyzed with the ChemStation software, and 
the relative abundance of different hydrocarbon residues was 
calculated by the ratio of the peak area of each hydrocarbon to the 
peak area of n-pentadecane in the GC chromatograph. The residue 
ratio of n-alkanes was calculated with the equation R—[S]/[I], 
where R, [S], and [i] represent the alkane residue ratio, the residual 
n-alkane concentrations in the samples, and the concentration before 
inoculation, respectively. The results were expressed as mean values 
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Table 1. Strains, plasmids and putative genes studied in the work. 





Name 


Description 


Source/reference 


Strains 


P. aeruginosa SJTD-1 


Host strain isolated from oil-contaminated soil, Wild Type 


This study 


Mut a i kB1 Mut a i kB2 Mut alkB1&2 


AlkB-deleted mutants of P. aeruginosa SJTD-1 


This study 


M u t P450 . 1 /Mutp4 50 . 2 /M Ut P450 . 1&2 


P450-deleted mutants of P. aeruginosa SJTD-1 


This study 


Mut a/mA 


AlmA-deleted mutants of P. aeruginosa SJTD-1 


This study 


Plasmids 


pRKaraRed 


Broad-host-range, X-Red proteins expression vector, Tet R 


[25] 


Genes' 


4712 


alkane-1 monooxygenase, AlkBI 


This study 


3623 


alkane-1 monooxygenase, AlkB2 


This study 


5482 


putative cytochrome P450 hydroxylase, P450-1 


This study 


4609 


putative cytochrome P450 hydroxylase, P450-2 


This study 


3206 


monooxygenase, flavin-binding family, AlmA-like protein 


This study 



*Gene numbers presented in RAST server. 
doi:1 0.1 371 /journal.pone.01 05506.W01 



with standard deviation. The cell-free controls were incubated and 
analyzed in the same manner. 

Deletion of the putative alkane-degradation genes from 
strain SJTD-1 

Based on the whole-genome DNA sequence of strain SJTD- 1 
[22], two putative alkB genes, gene 3623 and gene 4712, showed 
great similarities to the reported alkB2 and alkBl genes, 
respectively. To determine the function of alkBl and alkB2, the 
knockout mutants of the single gene and two genes were 
constructed using the pRKaraRed homologous recombination 
system [25]. Briefly, plasmid pRKaraRed was first electro- 
transformed into strain SJTD-1 to generate the SJTD- 1 /pRKar- 
aRed competent cells. The in-frame scarless knockout mutant of 
the single gene was generated by two-step homologous recombi- 
nation. First, the sacB-bla cassettes were amplified from plasmid 
pEX18Ap [26] with primer pairs Fl and Rl, and then electro- 
transformed into the SJTD-1 /pRKaraRed competent cells. The 
transformants were screened on LB plates supplemented with 
500 (Xg/mL carbenicillin and 50 Ug/mL tetracyclin. Next, the 
sacB-bla removal cassettes amplified from the genomic DNA of 
the first-step colonies (Suc s Carb R ) using the primer pairs F2 and 
R2, were electro-transformed into the competent cells of the first 
step to perform the second recombination. The transformants 
were spread on LB plates with 10% sucrose and 50 Ug/mL 
tetracycline. The transformants were further selected in parallel on 
LB plates with 10% sucrose and 50 (Xg/mL tetracycline, and on 
LB plates with 500 Ug/mL carbenicillin and 50 u.g/mL tetracy- 
cline. A positive genotype (Suc R Carb s ) was verified by PCR 
detection with the test primers and DNA sequencing was 
performed by Invitrogen. The deletion of multiple genes was 
achieved by repeating this two-step knockout method. Two alkB 
genes in strain SJTD-1 were knocked out separately or simulta- 
neously, and the mutant strains were named as Mut n a Bi , 
Mut a H£2, an d Mut,,MJ6?2, respectively. Similarly, two putative 
P450 genes (gene 5482 and 4609), and one AlmA-like monoox- 
ygenase gene (gene 3206) were also deleted and the mutant strains 
were Mutp4jo-l> ^^tp4so-2> an< J MutAimA, respectively. All the 
mutants are listed in Table 1; all the primers for homologous 
recombination and PCR detection are fisted in Table SI. The 



sequences of the five alkane hydroxylase genes are shown in Text 
SI. 

To study the roles of these alkane hydroxylases in n-alkane 
utilization, the growth curves of strain SJTD-1 and mutants (with 
different w-alkanes as the sole carbon sources) were determined by 
the Automatic Growth Curve Analyzer as described above. 
Experiments were replicated thrice and the results obtained were 
the averages values with standard errors. 

Quantitative reverse transcription-PCR (Q-RT-PCR) 

To study the transcription profiles of the five alkane hydrox- 
ylases in different conditions, the wild type and mutant SJTD-1 
were cultured in the BSM medium with sodium acetate and 
various alkanes (C 12 -C 2 4) to the mid-exponential phase. The 
alkanes were the sole carbon sources. Total RNA was extracted 
using the FastRNAPro Blue Kit for microbes (MP Biomedicals, 
Santa Ana, CA) and the RNA extraction was performed according 
to the manufacturer's instructions. The RNA yield was estimated 
using a Nanodrop UV spectrometer (Thermo Scientific, Wilming- 
ton, DE, USA). The reverse transcription was achieved using the 
PrimeScript Reverse Transcriptase Kit (TaKaRa, Dalian, China). 
Approximately 1 ug RNA and 20 ng random primers were used. 
The Q;RT-PCR was performed using the IQSYBR Green Super- 
mix (Bio-Rad Laboratories, Hercules, CA) and gene-specific 
primers (RT-F and RT-R, Table SI) in the IQTM 5 Multicolor 
Real-time PCR Detection System (Bio-Rad, CA, USA). The Q,- 
RT-PCR procedure carried out was as follows: 15 min of pre- 
denaturing, 40 cycles of 95°C for 10 s; 60°C for 30 s; and 72°C for 
30 s; followed by a melting curve stage from 60°C to 95°C. The 
cDNA amplification efficiency of the samples, internal standards 
(16S rRNA), and calibrators (samples without induction by 
alkanes) were equivalently modulated, and the relative fold change 
in mRNA quantity was calculated using the DDCt method [27]. 
At least three independent Q,-RT-PCR experiments were 
conducted for each RNA sample. The change in the transcrip- 
tional levels of strains grown with alkanes was calculated by 
comparing with the transcriptional level of the strains grown in 
sodium acetate [17]. 
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5G 



50 
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Pseudomonas chlororaphis subsp. chlororaphis DSM 50083T (NR_044974.1) 
— Pseudomonas chlororaphis subsp. aureofaciens DSM 6698 (NR_042939.1) 

Pseudomonas libanensis CIP 105460 (NR_024901.1) 

Pseudomonas fluorescens 1AM 12022 (NR_043420.1) 



Pseudomonas mandelii CIP 105273 (NR_024902.1) 

■Pseudomonas rhizosphaerae IH5 (NR_029063.1) 



50 



Pseudomonas putida 1AM 1236 (NR_043424.1) 

- Pseudomonas plecoglossicida FPC951 (NR_024662.1) 



100 I Pseudomonas putida GPol (AJ249825.2) 

Pseudomonas xanthomarina KMM 1447 (NR_04 1044.1) 



99 



- Pseudomonas fiavescens B62 (NR_025947.1) 

Pseudomonas argentinensis CHOI (NR_043115.1) 




90 



70 



Pseudomonas mendocina NCIB 10541 (NR_043421.1) 

Pseudomonas pseudoalcaligenes Stanier 63 (NR_037000.1) 
Pseudomonas alcaliphila AL15-21 (NR_024734.1) 

Pseudomonas alcaligenes IAM12411 (NR_043419.1) 

Pseudomonas oleovorans 1AM 1508 (NR_043423.1) 



100 



■ Pseudomonas otitidis MCC10330 (NRJ343289.1) 

Pseudomonas aeruginosa DSM 50071 (NR_026078.1) 



100 



99 



I Pseudomonas aeruginosa DQ8 (GU269267.1) 

- ▲ Pseudomonas sp. SJTD-1 (JQ951926.1) 



69 L Pseudomonas aeruginosa PAOl (DQ.777865.1) 



0.005 



Figure 1. Phylogenetic tree based on 16S rDNA gene sequences indicating the SJTD-1 position. The Kimura two-parameter distance 

model was used, and bootstrap analysis was performed with 1,000 repetitions using MEGA 5.0. 

doi:10.1371/journal.pone.0105506.g001 



Results 

Isolation and identification of the oil-degrading strain 
SJTD-1 

Strain SJTD-1 isolated from the enrichment culture of oil- 
contaminated soil is a rod-shaped, gram-negative bacterium. Its 
growth pH ranged from 4.0 to 10.0, and optimal growth occurred 
at pH 7.0. After growing on LB agar at 37°C for 24 h, the cells 
formed green, round, moist and glossy colonies, approximately 
1.0 mm in diameter. Strain SJTD-1 could utilize n-alkanes, from 
n-dodecane (C 12 ) to n-triacontane (C 30 ), as its (sole) carbon 
sources. Although SJTD-1 could not grow in the presence of n- 
decane or other shorter length n-alkanes, its growth was not 
inhibited by the short-chain n-alkanes like n-hexane. The 16S 
rRNA gene sequence of SJTD-1 (Genebank Accession No. 
JQ951926.1) was 99% identical to that of P. aeruginosa PAOl 
(Genebank Accession No. DQ777865.1), and the corresponding 
phylogenetic tree analysis supported a strong relationship between 
SJTD-1 and members of Pseudomonas sp. (Fig. 1). Therefore, 
SJTD- 1 was classified as a P. aeruginosa strain. 

Growth curve of strain SJTD-1 and its n-alkane 
degradation efficiency 

The SJTD-1 -mediated breakdown and utilization of n-alkanes 
was studied by monitoring its cell growth for seven days in a BSM 



medium containing a variety of n-alkanes (Ckj-C^) at various 
concentrations (100 mg/L to 2 g/L) (Fig. 2). Cells started 
multiplying immediately after incubation and reached the 
exponential growth phase between day 1 and day 3. Of all the 
n-alkanes, n-hexadecane caused the highest amount of cell 
growth, followed by n-octadecane and n-eicosane, thus suggesting 
that n-hexadecane was probably the best available carbon source 
for strain SJTD-1 in this set of experiments. At higher 
concentrations of n-alkanes, strain SJTD-1 accumulated more 
biomass (as reflected by OD 60 o), and also required more time to 
reach the stationary phase. The stationary phase was attained after 
36 h when 500 mg/L of n-tetradecane, n-hexadecane, n-octade- 
cane, and n-eicosane were used. At 1 g/L of n-dodecane and n- 
docosane, the stationary phase was reached after 60 h and 36 h, 
respectively. When 2 g/L n-tetracosane was used, the stationary 
phase was attained after nearly 100 h (Fig. 2). 

We also determined the n-alkane degradation efficiency of 
strain SJTD-1 with GC-MS. The corresponding time courses are 
illustrated in Fig. 3. GC-MS results showed that substrate 
consumption was closely followed by biomass increase (the loss 
due to vaporization and extraction was excluded by running two 
blank controls). As shown in Fig. 3A to 3C, more than 95% of 
500 mg/L n-tetrodecane, n-hexadecane, and n-octadecane were 
degraded within 36 h (from 500 mg/L to about 15 mg/L, 5 mg/ 
L, and 5 mg/L, respectively). Complete degradation occurred in 



PLOS ONE | www.plosone.org 



4 



August 2014 | Volume 9 | Issue 8 | e105506 



Alkane Degradation Strain SJTD-1 and Its Alkane Hydroxylase Genes 






PLOS ONE | www.plosone.org 



5 



August 2014 | Volume 9 | Issue 8 | e105506 



Alkane Degradation Strain SJTD-1 and Its Alkane Hydroxylase Genes 



Figure 2. Growth curves of P. aeruginosa SJTD-1 in /7-alkanes. Strain SJTD-1 was cultured in BSM supplemented with different concentrations 
(■ 100 mg/L, • 250 mg/L, ▲ 500 mg/L, T 1,000 mg/L, and ♦ 2,000 mg/L) of n-alkanes: n-decane (A), n-dodecane (B), n-tetradecane (C), n- 
hexadecane (D), n-octadecane (E), n-eicosane (F), n-docosane (G), and n-tetracosane (H) for seven days at 30°C. The data are expressed as mean 
values and the standard deviations are shown. 
doi:1 0.1 371 /journal.pone.01 05506.g002 



less than 48 h, and the biomass at OD 600 of n-tetrodecane, n- 
hexadecane, and n-octadecane increased to 0.35, 0.72, and 0.50, 
respectively. It was obvious that the degradation rates of strain 
SJTD-1 for various lengths of n-alkanes were different. Strain 
SJTD-1 could degrade n-eicosane, n-docosane, and n-tetracosane 
completely until the end of the experiments (96 h or 160 h). Only 
about 50% of the n-eicosane was degraded in 96 h. Although a 
rapid decline to about 70% residues was observed in the first 24 h, 
the degradation rate decreased in subsequent hours and the net 
biomass accumulation corresponded to an OD 6 ( )0 ~0.58 (Fig. 3D). 
Similarly, the utilization efficiencies for n-docosane and n- 
tetracosane were about 35% and 25%, respectively, and their 
respective biomass reached an OD 6m of approximately 0.39 and 
0.26, respectively (Fig. 3E and 3F). All these indicated that strain 
SJTD-1 degraded the medium-chain n-alkanes more rapidly than 
the long-chain n-alkanes, probably due to the much higher 
hydrophobicity of the latter [17,28,29]. 

Moreover, to find the highest concentration of n-alkanes that 
the SJTD-1 strain could tolerate and utilize for growth, we 
analyzed its degradation efficiency for n-hexadecane (concentra- 
tion ranging from 250 mg/L to 2.0 g/L) because n-hexadecane 
was possibly its greatest carbon source. As shown in Fig. 4, all the 
n-hexadecane could be completely degraded in three days. 
Additionally, 250 mg/L and 500 mg/L of n-hexadecane were 
totally degraded within 36 h; as for 1.0 g/L and 2.0 g/L n- 
hexadecane, approximately one more day was needed for a 
complete degradation. Therefore, we concluded that SJTD-1 
could tolerate and completely biotransform 2.0 g/L n-hexadecane 
into a relatively large amount of biomass (OD 600 = 1.7, Fig. 2). 
This result reveals its great tolerance ability, fast degradation 
speed, and high utilization efficiency for n-alkanes. 

Sequence analysis and alignment of the alkane- 
degradation genes in strain SJTD-1 

The whole-genome DNA sequence of strain SJTD- 1 has been 
obtained and deposited in GenBank under the accession number 
AKCM00000000 [23]. It was also annotated by the RAST server. 
Two putative genes encoding AlkB monooxygenases were found 
to be involved in the first step of alkane degradation. Gene 47 1 2 
(alkB 1) encodes a putative AlkBl monooxygenase, comprising 383 
amino acids, and with an estimated molecular weight of 42.1 kDa. 
The second alkB gene 3623 (alkBT) encodes AlkB2 monooxygen- 
ase, comprising 378 amino acids and with an estimated mass of 
41.6 kDa. The sequence oialkBl and alkBI genes exhibited 78% 
identity, and their derived proteins exhibited 68% identity. The 
cross-genera sequence alignments of AlkBl and AlkB2 with other 
published AlkB sequences indicated that all of these alkane 
hydroxylases share one conserved HYG motif and three highly 
conserved regions of the Hist boxes containing eight histidines 
[30-32] (Fig. 5A). 

To further understand the relationship between AlkB proteins, 
phylogenetic analysis was performed based on the characterized 
AlkB monooxygenases from Pseudomonas, Alcanivorax, Acineto- 
bacter, Mycobacteria, and Rhodococci (Fig. 5B). In the phylogenetic 
tree, the two AlkB proteins from P. aeruginosa PAOl and Ml 8 
were clustered as alkane monooxygenase 1 , and these were close to 
the alkane monooxygenases from other gram-negative bacteria, 



such as Alcanivorax borkumensis SK2 and Acinetobacter sp M-l. 
In fact, the AlkBs of strain SJTD-1 showed very high sequence 
identity with other Pseudomonas AlkBs. For instance, the identity 
of SJTD-1 AlkBl to P. aeruginosa PAOl and Ml 8 AlkBl was 
100%, and the AlkB 2 identity was 100% and 99%, respectively. 
Another reported Pseudomonas strain DQ8, that could easily 
break down alkanes and polycyclic aromatic hydrocarbons 
harbored only one AlkB homologue, which was mapped to 
SJTD-1 AlkB 2 with 100% identity [8]. In addition, two 
cytochrome P450 monooxygenase homologues (gene 5482 and 
4609) and one AlmA-like monooxygenase homologue (gene 3206) 
were also found in the genome of strain SJTD-1. Sequence 
analysis showed that genes P450-1 (gene 5482) and P450-2 (gene 
4609) encoded a 418-amino-acid-protein and a 444-amino-acid- 
protein, respectively; they shared 99% amino acid sequence 
identity to those of P. aeruginosa PAOl and 26% to those of 
Alcanivorax dieselolei B-5 [33]. As for the putative almA-like gene, 
its protein (499 aa) displayed 56% identity with the putative flavin- 
containing monooxygenase of A. dieselolei B-5, and only 50% 
identity to the AlmA of Acinetobacter sp. DSM 17874 [34]. 

As shown in Fig. 5C, the genomic organization of the five 
alkane hydoxylases, none of which are located within the gene 
cluster encoding the genes responsible for alkane catabolism like 
those seen in P. putida GPol and other hydrocarbonoclastic 
bacteria [19,33,35]. It should be noted that the surrounding open 
reading frames of the alkB genes in strain SJTD-1 were different 
from those in other strains, while the rubredoxin (Rd)-encoding 
genes, the rubredoxin reductase-encoding genes, and the tran- 
scriptional regulatory protein-encoding genes were often located 
immediately downstream of the alkB gene. No rubredoxin 
reductase-encoding gene was sought out; and no fusion Rd 
domain was found near the alkB genes, as opposed to that 
observed in Dietzia sp. DQ12-45-lb [31]. However, some 
important clues were found despite the decentralized distributions 
over the chromosome [36]. One exception was tlpS, a gene coding 
for a methyl-accepting chemotaxis protein that is potentially 
involved in chemotaxis towards long-chain n-alkanes. It was found 
located downstream of alkB2 [37]. A putative transcriptional 
regulator of the GntR family was located upstream of AlkB2 and 
the intergenic region between alkB2 and gntR was only 181 bp. 
Therefore, it is conceivable that in strain SJTD-1, n-alkanes may 
be utilized through a more diversified pattern than that in other 
reported hydrocarbonoclastic bacteria. 

Phenotypic study of the putative alkane hydroxylation 
genes in strain SJTD-1 

In order to determine the functions of the two alkB genes in the 
n-alkane oxidative process and in order to understand their 
preferences to n-alkanes, we knocked out the two genes separately 
and simultaneously, and then monitored the utilization capability 
of the three mutants (Mut a a#i, Mut a ffi#2, and Mutaft^;^) for the 
C 12 -C 2 4 n-alkanes. The mutants were verified using PCR 
detection (Fig. SI). For the medium-chain n-alkanes (C 12 -C 16 ), 
single-knockout mutants showed an obvious delay as compared to 
the wild type SJTD-1, and double-knockout mutants could not 
utilize the n-alkanes at all (Fig. 6A-6C). It meant that these two 
genes were responsible for the degradation of C 12 -C 16 n-alkanes 
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Figure 3. Growth of strain SJTD-1 in different n-alkanes (500 mg/L) and its consumption of n-alkanes. ■ indicated as OD 600 in the left 
vertical axis; Aexpressed as n-alkane residues detected by GC-MS with respect to abiotic controls in the right vertical axis); A indicated as the n-alkane 
residues of the abiotic controls. A, B, C, D, E, and F represented the n-alkanes C n4 , C q6 , C q8 , C 2 o, C 2 2, and C 2 4, respectively. Standard errors were 
calculated from three independent determinations. 
doi:10.1371/journal.pone.0105506.g003 



and that both of them were involved in the oxidative process. The 
alkB2 defect showed a more pronounced effect on the cell growth 
than did the alkBl defect, thereby suggesting that AlkB2 played a 
more major role in the oxidation of medium-chain alkanes 
(Fig. 6A-6C and Fig. 7). Unexpectedly, the poor viability of 
Mut„ftBi&f2 recovered to a normal state when C 18 -C 2 4 alkanes 
were used as the sole source of carbon. No obvious growth 
decrease was found in the single-knockout mutants (Fig. 6D-6G). 
This implied that there were some hitherto unknown genes 
involved in the utilization of longer-chain n-alkanes. 

We further studied three more potential alkane hydroxylases 
(AH) genes (two putative cytochrome P450 genes and one almA 



monooxygenase gene) listed in Table 1. In order to judge the 
physiological role of each putative AH, single-knockout, double- 
knockout or multi-knockout mutants were obtained, and their 
viabilities were evaluated in the BSM media supplemented with n- 
alkanes (Fig. 7). On the whole, the divergence of the mutants' 
viabilities in Fig. 7 resulted from the dissimilar enzymatic 
preferences of AHs within the range of the substrate spectrum. 
Strains with absence of two AlkBs almost failed to grow on C^— 
C 16 alkanes; whereas P450 and AlmA made a minor contribution 
to the degradation of these substrates, as their defective strains still 
grew to a comparable level relative to the parent. However, it 
could be concluded that the recovery of viability originated from 
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Figure 4. Quantitative estimation of the SJTD-1 -mediated degradation of /r-hexadecane. The strain was cultured in BSM supplemented 
with 250 mg/L (■), 500 mg/L (•), 1,000 mg/L (A), and 2,000 mg/L (T) of n-hexadecane at 30 r C for seven days; A indicated as the n-alkane residues 
of the abiotic controls. Standard errors were calculated from three independent determinations. 
doi:10.1371/journal.pone.0105506.g004 



the transcription of P450 and AlmA oxygenases, as mutant 
M-Ut n // !B1 tf2-P450-i&2-almA showed a similar growth level as that of 
the wild type for C lf .-C 2 4. Therefore, an affirmative hypothesis for 
the contribution of other undefined AHs to microbial catabolism 
of long-chain alkanes (Cig— C24) was made. So far, these findings 
explain the metabolic roles of identified AHs within their preferred 
substrate range and further provided clues for exploring the 
uncharacterized alkane hydroxylases. 

Transcriptional expression profiles of the alkBI, alkB2, 
p450 and almA genes in strain SJTD-1 

Five putative AHs in strain SJTD-1 were revealed as expected, 
and the diversity of AHs allowed for a rapid adaptation of SJTD-1 
to the changing environments of various alkanes. To articulate the 
contribution of each AH to the utilization of alkanes, their 
transcriptional profiles were investigated. Fig. 8A shows that C12— 
C 2 4 alkanes rendering the viability of SJTD-1 rested with a 
universal expression of AHs in any case. AlkBI showed a relatively 
narrow inducer range of C12-C16 than did AlkB2, which could be 
induced by Ci 2 -C 2 4 alkanes, thereby eliciting more than a two- 
fold change on its transcriptional level (Fig. 8A). Two P450 
monooxygenases probably contributed to the mineralization of 
C12-C16 alkanes, because of the 2.5-fold change with those substrates. 
However, when Ci 8 -C 24 alkanes were used, their transcriptional 
changes fell down sharply (Fig. 8A). The AlmA-like oxidase was 
expressed at a low level (1.5-3 fold change), only when the Ci 8 -C 2 4 
n-alkanes were used (Fig. 8A). So far, it was facile to explain that the 
abundant transcription of AHs in the hexadecane-induced cells lead 



to the most significant utilization of n-hexadecane during SJTD- 1 
growth in Fig. 2. 

As the single deletion of the five AHs made little impact on the 
cell growth, the transcription levels of the targeted AHs were 
analyzed in the alkB 1 -defective, aftB2-defective, and double- 
deletion mutants. A functional complementation of AHs was 
observed during the cultivation period. As illustrated in Fig. 8B 
and 8C, the transcriptional level of an alternative AlkB hydrox- 
ylase was enhanced and it was thus responsible to terminally 
oxidize n-alkanes, particularly when one AlkB was absent. This 
terminal oxidation of w-alkanes was limited as indicated by the 
overlapping area of both AlkB substrate spectrums, i.e. C 12 -C 16 
alkanes. Among these substrates, two P450s also displayed changes 
in expression and contributed to the complementary effects of 
AlkBs. Interestingly, the AlmA-like oxygenase increased its 
transcription to some extent (within the range of C ls -C 24 alkanes) 
when alkB2 was removed rather than alkBI (Fig. 8B-8D). 
Therefore, the less attenuated viability of strain SJTD-1 on 
alkanes even when any alkB gene was inactivated can be attributed 
to the upregulation of functionally-related AHs and their 
metabolic contribution. 

As mutants lacking both the AlkB hydroxylases showed poor 
propagation on C12-C16 alkanes, but did not affect the duplication 
rate on the C 18 alkane, the transcription of two P450 and AlmA- 
like oxygenases was surveyed as depicted in Fig. 8D. However, no 
obvious alteration was found for these AHs in the single or double 
alkB mutants. Mutants with a consecutive deletion of five AHs still 
survived on n-octadecane and also duplicated to a normal level of 
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Note: sthk: Signal transduction histidine kinase; traits: Transporter; 
fadA: Acyltransferase; xdh-FeS: Xanthine dehydrogenase iron-sulfur; 
xdh-MoII: Xanthine dehydrogenase Mo(II) binding; zipA: Cell division protein; 
Ars:3-Oxoacyl-ACP synthase; araJ: Arabinose efflux permease; 
ftr:3-Oxoacyl-ACP reductase; zetA: Glutathione s-transferase; 
dsbD,E,G: Cytochrome c-type biogenesis protein; 
kefB: Glutathione regulated potassium efflux system protein; 



Figure 5. Sequence alignment of alkane hydroxylases AlkEtl and AlkB2 from several strains and genetic distribution of five putative 
alkane hydroxylases in P. aeruginosa SJTD-1. Strains list: SJTD-1, P. aeruginosa SJTD-1; PA01, P. aeruginosa PAOI; DQ8: P. aeruginosa DQ8; GPol, P. 
putida Gpol; PI, P. putida PI; CHAO: P. fluorescens CHAO; SK2, Alcanivorax borkumensis SK2; MA , Acinetobacter sp. M-1; Q15, Rohodococcus erythropolis 
Q1 5; H37Rv, Mycobacterium tuberculosis H37Rv; DQ1 2-45-1 b, Dietzia sp. DQ1 2-45-1 b; E1 , Dietzia sp. E1 . (A) The conserved amino acid residues in all the 
sequences are indicated with a blue background, and the dark degree is in direct proportion to the degree of conservation. The three conserved 
histidine boxes (Hist-1, Hist-2, and Hist-3), and one HYG motif are underlined. The conserved amino acid associated with each position is shown below 
the alignment (created by Clustal Origin). (B) Phylogenetic relationship based on the complete amino acid sequences of AlkB1 and AlkB2 from P. 
aeruginosa SJTD-1 and other published AlkB amino acid sequences. The phylogenetic tree was determined using the Kimura two-parameter distance 
model, and bootstrap analysis was performed with 1,000 repetitions using MEGA 5.0. (C) Schematic profile of the genomic organization of AHs 
homologues in P. aeruginosa SJTD-1 . alkBl and alkB2, alkane-1 monooxygenases (lawngreen and grey, separately); P450-1 and P450-2, cytochrome P450s 
(red and blue, respectively); almA, AlmA-like monooxygenase (green); tlpS, methyl-accepting chemotaxis protein. 
doi:1 0.1 371 /journal.pone.01 05506.g005 
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Figure 6. Growth curves of P. aeruginasa SJTD-1 (■), Mut, 1/kB , (•), N\ut a i kB2 (A), and Wlut a / tB7 ^ (T) in n-alkanes. The strains were inoculated 
in BSM supplemented with 500 mg/L n-alkanes: n-dodecane (A), n-tetradecane (B), n-hexadecane (C), n-octadecane (D), n-eicosane (E), n-docosane 
(F), and n-tetracosane (G) and cultured for seven days at 30°C. 
doi:10.1371/journal.pone.0105506.g006 
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Figure 7. 3D-bar growth representations of all the defective mutants with C 12 -C 24 n-alkanes. Values on the horizontal axis profile the 
genotype of the mutants with single or multiple deletions of alkane hydroxylases; The Y-axis specifies the tested n-alkanes, coloring from black to 
dark blue black. The Z-axis of bar altitudes represents the OD 60 o for the viability of each mutant. Sharp contrasts were symbolized using arrows when 
both alkBl and alkB2 were knocked out, and cultivated in C 16 and C q8 alkane media, respectively. 
doi:1 0.1 371 /journal.pone.01 05506.g007 



ODgoo (Fig. 7). This strongly suggested that the biodegradation of 
Cqg— C 2 4 alkanes differed from that of the C 12 -C 16 medium chain 
alkanes, and that a cluster of unidentified AHs may participate in 
this process, by transforming the long chain alkanes into lower 
hydrocarbons. 

Discussion 

P. aeruginosa strain SJTD-1, a new soil isolate, was determined 
as a hydrocarbonoclastic consumer of medium- and long-chain 
alkanes (G12-C30). Based on the 16S rDNA sequence alignments 
and phenotype analysis, strain SJTD-1 showed close homogeneity 
to the well-documented P. aeruginosa PAOl and P. aeruginosa 
DQ8. However, strain SJTD-1 exhibited no human pathogenicity 
like PAO 1 , and it degraded polycyclic aromatic hydrocarbons only 
to some extent, unlike DQ8. In fact, a relatively narrower substrate 
range of C 12 -C 24 supported a sustainable growth of strain SJTD- 
1. For this range of n-alkanes, strain SJTD-1 showed outstanding 
degradation efficiency and utilization speed. n-Hexadecane 
followed by n-tetradecane and n-octadecane were assimilated as 
preferential carbon sources for this strain as these hydrocarbons 
could be rapidly biodegraded by SJTD-1. M-Tetradecane, n- 
hexadecane, and n-octadecane (at a concentration of 500 mg/L) 



were completely degraded within 36 h. SJTD- 1 could also tolerate 
and degrade 2.0 g/L M-hexadecane within three days. As 
compared to other alkane-utilizing strains, strain SJTD-1 degrades 
the same type and same amount of n-alkanes more efficiently and 
more quickly. Most of the alkane-consuming strains mineralize 
equivalent amounts of alkanes at a much slower rate and require a 
much longer degradation time (Table 2). When 500 mg/L of n- 
octadecane was depleted by strain SJTD-1 in 36 h, Alcanivorax 
sp. 2B5 simply consumed 16.07% of this substrate at the beginning 
of 48 h. In contrast to P. aeruginosa DQ8, strain SJTD-1 showed 
a comparable transformation velocity but the initial concentration 
of n-tetradecane was 5-fold higher in this case. 

In order to decipher the mechanism by which SJTD-1 breaks 
down the n-alkanes, the initial but rate-limiting step deserves 
special attention. This step marks the beginning of the terminal 
oxidation catalyzed by alkane hydroxylases. The genome sequenc- 
ing results indicated that at least five alkane hydroxylases in strain 
SJTD-1 were involved in the oxidization process. Through 
transcriptional analysis, two alhB genes were found to be gready 
induced by n-alkanes, with AlkB2 sharing a broader substrate 
range. The two P450 monooxygenases were induced only when 
medium-length alkanes were used, whereas the Ci 8 -C 2 4 alkanes 



PLOS ONE I www.plosone.org 



11 



August 2014 | Volume 9 | Issue 8 | e105506 



Alkane Degradation Strain SJTD-1 and Its Alkane Hydroxylase Genes 




NaAc C12 C14 C16 C18 C20 C22 C24 NaAc C12 C14 C16 C18 C20 C22 C24 NaAc C12 C14 C16 C18 C20 C22 C24 



Figure 8. Transcriptional levels of five putative AHs genes in wild type SJTD-1 and in three knockouts. Relative expression levels were 
determined with real-time RT-PCR. 16S rRNA gene was used as the reference gene. The expression levels of genes in strain SJTD-1 (A), strain Mutate, 
(B), strain Mut ottB2 (C), and strain Mut a i kB , &2 (D) were shown. All strains were grown in different n-alkanes (C 12 to C 24 alkanes as the sole carbon 
source), and cells grown with sodium acetate were used as controls. Standard errors were calculated from three independent determinations. 
doi:10.1371/journal.pone.0105506.g008 



were able to regulate the expression of the AlmA oxygenase, which 
is known to be involved in the oxidation of super long-chain 
alkanes (>C 30 ) [34]. The fact that w-hexadecane showed a 
prominent growing effect for strain SJTD-1 was evident from the 
topmost transcription of AHs during hexadecane-induced culture, 
as compared to the other n-alkanes. A transcriptional enhance- 
ment of two AlkBs was evident when either of the AlkBs was 
inactivated, strongly suggesting their overlapping substrate spec- 
trum and functional identities. 

In order to specify the roles of AHs in the alkane degradation, 
some important defective strains were deliberately constructed and 
their phenotypic profiles were represented under the given range 
of C 12 -C 24 alkanes. For medium chain alkanes (C 12 -C 16 ), AlkB 
hydroxylases showed a dramatic performance, and bacterial 
growth was seriously influenced when the two alkB genes were 
deleted. Moreover, the defective strain in absence of AlkB2 grew at 
a slower rate than the AlkB 1 defective mutant did. It was therefore 
speculated that AlkB2 played a central role in the catabolism of 

Table 2. Comparison of alkane degradation by different strains. 



C 12 — Ci 6 alkanes, which was more important than the role played 
by AlkBl. The other three AHs including P450 and AlmA 
oxygenases had little effect on n-alkane degradation, although 
their transcription would somewhat offset the deletion of the alkB 
genes. This implies a heavy redundancy of AHs for the microbial 
degradation of w-alkanes. However when Ci 8 -C 24 n- alkanes were 
used, strain SJTD-1 utilized another set of alkane hydroxylases to 
metabolize these n-alkanes. Therefore, it is important to under- 
stand how strain SJTD- 1 , in the absence of all known AHs, still 
manages to survive on the long-chain alkanes. 

The key factors involved in the degradation of alkanes, such as 
AlkB, rubredoxin, and other related alkane hydroxylases, showed 
a scattered distribution in the genome of strain SJTD-1, without 
the occurrence of gene clusters as represented in Pseudomonas 
putida cells. This was also observed among other P. aeruginosa, 
such as P. aeruginosa PAOl and P. aeruginosa RR1. However, 
their sporadic distribution did not compromise the bacterial 
utilization of n-alkanes. This indicates that a compensatory 





Strain 


Substrates 


Conservation rate 


Sources/Reference 


P. aeruginosa SJTD-1 


Cl4< ^-16/ C 18 


>99% in 36 h 


This study 


500 mg/mL 


P. aeruginosa DQ8 


C, 4 , 100 mg/L 


>99% in 24 h* 


[8] 


Alcanivorax sp. 2B5 


C, 8 , 500 mg/L 


1 6.07% in 48 h 


[32] 


>99% in 96 h 


Acinetobacter 


Cie, 400 mg/L 


>99% in 60 h 


[39] 


Dietzia DQ1 2-45-1 b 


C 16 , 2.5 mg/mL 


31% in 21 d 


[11] 



*with 0.005% yeast extract in the medium. 
doi:1 0.1 371 /joumal.pone.01 05506.t002 
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mechanism was responsible for the regulation and coordination of 
these multi-step catalytic reactions. For instance, it has been 
proposed that genes encoding two significant electron transfer 
proteins, viz. rubredoxin and rubredoxin reductase, were distant 
from the inducible alkB genes and constitutively expressed in any 
case [38]. To better understand the mechanism behind this 
disordered genetic organization, we need to resort to a compre- 
hensive approach, such as an iTRAQ-based proteomics study, 
which would systematically characterize the global response of 
strain SJTD-1 to various types of medium- and long-chain n- 
alkanes. 

Supporting Information 

Figure SI PCR verification of the mutant strains. The 

genomic DNA of all the mutants as well as the wild type P. 
aeruginosa SJTD-1 was extracted and amplified with the test 
primers listed in Table SI. The colonies with shorter PCR 
products (400-500 bp) represented successful deletions. (A) The 
PCR products of Mut aft£i , Mut a ikBi&2> and wild type SJTD-1; (B) 
The PCR products of MuW™^ and wild type SJTD-1; (C)The 
PCR products of Mutp45o_;, Mutp4jo-2, Mutp45o ;&f2, and wild 
type SJTD-1. 
(TIF) 
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